o 

(N 

< 

(N 
(N 



O 



Stability and convergence of a higher order rational 

difference Equation 

Hamid Gazor*, Saeed Parvandeh"^ 



Department of Mathematics, Majlesi Branch, Islamic Azad University, Isfahan, Iran. 

Abstract 



Q I In this paper the asymptotic stabihty of equihbria and periodic points of the 

(-H ■ following higher order rational difference Equation 

» ' 

Xn+1 = -T— ' k>l, n = 0, 1, . . . 

p + ^^n^n—l ■ ■ ■ ^n—k 

is studied where the parameters a, /3, 7 are positive real numbers, and the initial 
^ ] conditions X-k, . . . ,xq are given arbitrary real numbers. The forbidden set of this 

equation is found and then, the order reduction method is used to facilitate the 
analysis of its asymptotic dynamics. 



Tj- , Keywords: Difference Equation; Equilibrium point; periodic point; conver- 

' gence; semiconjugacy 

o ■ 

1 Introduction and preliminaries 

^ ! Difference equations may appear as solutions of various phenomena or as a discretized 

I system of delay or non-delay differential equations. They are very important in both 

theory and applications; for applications in biology (see [2]), economics (see 0, E]), 
medical sciences (see [9]), military sciences (see The main goal in the study of 

difference equations is to understand the asymptotic behavior of solutions rather than 
trying to find an explicit formula for solutions. This basically is not only because the 
explicit solutions are hard to find but also the explicit solutions may still represent a 
complex dynamics and they yet may require a qualitative analysis to understand their 
dynamics. 

The ratio of any two polynomials of a recursive sequence is called a rational dif- 
ference equation. Rational difference equations are important as practical classes of 
difference equations. Most of the works in the literature of rational difference equations 
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have treated the first and second order difference equations. For second order rational 
difference equations we refer the reader to the monograph of Kullenovic and Ladas ([7]). 
In this paper we study the dynamics of the following {k + l)-order rational difference 
equation 

Xn+i = -5— , n = 0, 1, . . . (1) 

p -|- 'yXnXn—l ■ ■ ■ Xfi—k 

where k >1 and the parameters a,/3,7 are positive. We allow the initial conditions 
. . .xq to take any arbitrary value out of the forbidden set of equation ([1]). 
Let / be some interval of positive real numbers and / : I^^^ — )• J be a continuously 
different iable function. Then for every initial conditions (a;_fe, x_k+i) ■ ■ ■ , Xq) G 7*^"'"^, the 
difference equation 

Xn+l fi^-^nj ■^n—ly • • • ; ■^n— fc); 0, 1, . . . (2) 

has a unique solution The point x G / is called an equilibrium of equation ([2]) 

(or simply an equilibrium of /) if f{x for all n > (such a 

solution is also called a trivial solution). The point (ci, . . . , Cfc+i) is called a (A; + l)-cycle 
if = Cfc+i-j for alH = 0, 1, . . . , /c. In this case we say that {x„} is periodic 

with period {k + 1). The linearized equation associated with equation ^ about the 
equilibrium x is 

^n+l / ^ Q \X^ . . . XjZn—i, 
i=0 * 

and its corresponding characteristic equation is defined by 

An equilibrium x of equation(l2]) is called locally stable if for every e > 0, there exists 
6 > such that for the solution of equation (|2]) with \x-k—x\ + \x^k+i—x\+- ■ ■ + 

\xq—x\ < S we have x| < e for all n > 1. Furthermore, if there exists 7 > such that 
for the solution of equation ([2]) with — x| + x| + . . . + |xo— x| < 7 we 

have lim„_j.oo Xn = x, then x is called locally asymptotically stable. Linearized stability 
theorem indicates that if all roots of equation (jS]) lie inside the open unit disk |A| < 1, 
then the equilibrium x of equation ([2]) is locally asymptotically stable. If at least one 
of the roots of equation ([3]) has modulus greater than one, then the equilibrium x of 
equation ([2]) is unstable, see e.g., |6]. 

X is called a global attractor on an interval I if for every solution {xn}'^=-k of 
equation ([2]) with X-k,X-k+i, ■ ■ ■ ,xo G /, we have lim„^oo When X is locally 

stable and a global attractor, we call it globally asymptotically stable. The forbidden 
set of equation ([2]) is the set of all (k+l)-tuples (/i, /2, • • • fk+i), where they can not be 
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taken as the initial conditions for an infinite well-defined sequence {xn} on its domain. 
In other words, the forbidden set of equation ([2]) is the set of all initial conditions such 
that arbitrary iterations of the right hand side of equation ([2]) are not well-defined. 

The following definitions, lemma, corollary, and theorems are needed to study the 
global behavior of solutions of equation ([1]). 

Definition 1.1. (for original ideas see [3l [11]) Consider equation (|2]) and assume that 
D C ]R'^+^ is nonempty. 

(i) Suppose that (xi, ... 0:^+1) £ M'^^^ Then ||a;|| = max{a;i, . . . , x^+i} denotes the 

sup-norm of x. Also, if x is the equilibrium of / then X = (x, . . . ,x) G R'^"*'^ is an 
equilibrium of its vectorization that is 

Vf{xi,... Xk+l) = (/(Xi, . . . Xk+l),Xi, . . . , Xfc). 

(ii) If there is a non-constant function H G C(D,]R) such that HoVf = (j)oHonD 
for some G C(i/(D), M) then Vj is also called a [D, H, (fy-semiconjugate of M'^^-'^. 
The mapping H is called a link map and is called the factor map. For each 
t G H{D), the level set H~^{t) fl D, abbreviated as H^^ is called a fiber of H in 
D. 

(Hi) A continuous mapping /i : D — M is said to be bending at a point x G -D if x is 
not an isolated point in D and x ^ [/i^^(/i(x))]° fl-D, where 5"° denotes the interior 
of S. 

(iv) An equilibrium x of equation ([2]) is exponentially stable under / relative to some 
nontrivial interval / containing x if there is 7 G (0, 1) such that for every solution 
{x„} of equation ([2]) with initial values X-k, . . . , xq G / we have for all n > 1 that 
|x„ — x| < 07" where c = c{xq, . . . X-k) is independent of n. 

(v) An equilibrium x of a map / : M — t- M is semistable (from the right) if for any 

e > there exists 5 > such that if < Xq — x < 5 then |/"'(xo) — x| < e for all 
n > 1. If in addition, lim„_j.oo /""(^^o) = ^ whenever < Xq — x < 7 for some 7 > 0, 
then X is said to be semiasymptotically stable (from the right). Semistability 
(semiasymptotic stability) from the left is defined analogously. 

Theorem 1.1. (see 111])- Let Vf is a {D, H,(j))-semiconjugate map and x E D is an 
equilibrium ofVf. 

(i) t = H{x) is an equilibrium of (p. 

(ii) (Boundedness). Assume that |-f^(x)| —t- 00 as ||x|| — )■ 00. If the sequence {^'^(^0)} 
is bounded for some to £ H{D), then the sequence {Vp{xo)} with Xq G H^^ is 
bounded. 

(Hi) (Stability and instability) . Assume that H is bending at x. If we set t = H{x) 
then 
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(a) If X is stable (asymptotically stable) under Vj, then t is a stable (asymptoti- 
cally stable) equilibrium of (p. 

(h) Ift is unstable under cj), then x is unstable under Vf. 

(iv) (Attractivity of invariant fibers). Lett & I be an isolated equilibrium of (p which 
attracts all points in I. If Xq G Dr\H~^{I) and {Vp{xQ)} is bounded, then {Vp{xo)} 
converges to the invariant fiber (J). 

Theorem 1.2. (see IS]). Let x be an equilibrium of f : ^ M., fix) = 1, and 
f"(x) 7^ 0. Then x is semiasymptotically stable from the right (left) if f"(x) < 
(fix) > 0). 



2 The Forbidden Set 

Consider equation ([T]). If a = 0, the solution is trivial. If 7 = then equation ([T]) reduces 
to a linear equation. If /3 = then the solution is periodic with period k + 1. Thus, 

we assume that all the parameters are nonzero. A change of variables Xn = '^^\J\yn 
followed by the change yn = x„ reduces equation ([1]) to 

Xn+i = T— ^^^^^ , n = 0, 1, . . . (4) 

i + XfiXn—l ■ ■ ■ Xfi—k 

where c = ^ Hence, , we consider equation (jlj) instead of equation ([T]), hereafter. 

Now we are ready to obtain the forbidden set of equation (jlj). 

Theorem 2.1. Consider equation ([^. Assume that ^ is the forbidden set for this 
equation. Then 



00 

^ = [J I (/l; /2, • • • , /fc+l) : /1/2 • • • 



-1 



e: 



m=0 ^j=0 

Proof. Assume the initial conditions X-k,X-k+i, . . .xq satisfy 

-1 

for some m G N U {0}. Then equation (jlj) implies that 

-1 



_ cxqX^i . ..Xi_kX-k _ VErio^', 

XiXq . . . Xi-k — 1 I ^ „ ~ 1 v^m , — v-^m-l 



i + XQ...x-k 1 - YJIU 



j=0 



Therefore, we obtain from equation (jlj) that 



cxiXQ...X2-kXi-k _ '^VE™o''='y' _ ~1 



X2X1 . . . X2-k - , , - . . - „^_2 
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Continuing in this fashion by induction we obtain that 

_ -1 _ 

Ei=o c 

As a result the iteration process stops at x^+i- Conversely, assume iteration process 
stops at some point Xn^+i, i.e., Xn^Xno-i ■ ■ ■ Xn^-k = — 1. Then equation (jl]) implies that 

-1 

•^no — I'^riQ— 2 • • • ^riQ— A;— 1 | • 

Again using this fact and equation ^ we obtain that 

-1 

•^no— 2.^710— 3 . . . •^nQ—k—2 



1 + C + C2 ■ 

Continuing in this manner we get 

-1 

XqX-i . . .X-k 



or equivalently . . .Xq) G This completes our inductive proof. □ 



3 Linearized Stability 

The first step to understand of the dynamics is to find the equilibria and find their 
stability type. In this section we investigate the local asymptotic stability of the equi- 
libria of equation (jl]). Simple calculations show that origin is always an equilibrium for 
equation (jlj) and if c > 1 then it has a second equilibrium x = ''^\/ c — 1 . The following 
theorem deals with the local asymptotic stability of the equilibria. Note that throughout 
the rest of this paper, the initial conditions are assumed to be taken out of the forbidden 
set ^. 

Theorem 3.1. Consider equation ([^. Then, 

(a) origin is locally asymptotically stable for c < 1, while it is unstable when c > 1. 

(b) For the case of c> 1, the positive equilibrium x = '''^y/ c — 1 is stable. 

Proof. For part (a), the linearized equation associated with equation (jl]) about origin 
is Zn+i — cZn-k = 0, u > 0. Therefore, the corresponding characteristic equation is 
Xk+i _ c = 0, i.e., A = '=Vc. Since the origin is locally asymptotically stable if c < 1 
and is unstable if c > 1. 

(b) Simple calculations show that the linearized equation associated with equation (jl]) 
about the positive equihbrium x = '''y/ c — 1 is 

c- 1 c- 1 c- 1 1 _ _ 

Zn+l H Zn H Zn-1 + . . . H Z„_fc+i ^n-A: — U, n — U, 1, . . . 

CO CO 
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Therefore the corresponding characteristic equation is 

c-l c - 1 



A'=+i + A'^ + . . . + A - - = 0, (5) 

c c c 



Some algebra show that equation ([5]) is equivalent to 

(A'=+i-l)(A-l/c) 



A - 1 



0, A ^ 1, (6) 



Therefore, A = 1/c is one of the roots of equation ([6]). Since c < 1, by local asymptotic 
stability Theorem we conclude that the positive equilibrium x = V c — 1 is unstable. 
The proof is complete. □ 

Remark 3.1. The above theorem is sufficient to completely describe the local dynamics 
of the equilibrium for the case c < 1. However, for the case c < 1, more investigation 
are needed. For c = 1 the characteristic equation about origin has modulus equal to one, 
the linearization fails to analyze the asymptotic stability of origin. Thus, it may have a 
complex dynamics near this point. On the other hand consider equation ^ and assume 
that c > 1. equation ^ has a real root A = 1/c with modulus less than one but this 
equation has k + 1 roots in the following form 



2m7T \ ( 2m7T \ 

+ « sm 

k + lj \k + lj 



Xm = COS ( + i sin ( , m = 0, 1, . . . , /c, 



all with modulus equal to one. So, if c > 1 then linearization tells us nothing about the 
stability of the positive equilibrium. In the next section we discuss these cases in details. 



4 Semiconjugate factorization 

The main purpose of this section is to analyze the local dynamics near equilibria for the 
cases c > 1 as well as of the global asymptotic dynamics of the equation 01]). Although 
the previous section, using linearization, showed that the origin is locally asymptotically 
stable for c < 1. Yet the global nature of asymptotic stability can not be inferred from 
linearization. Also, we saw that for c > 1 the linearized equation about the positive 
equilibrium has several roots; all with modulus equal to unity. As a result, more powerful 
methods are necessary in order to analyze the dynamics of equation (jl]). In this section, 
we apply semiconjugacy analysis to examine the global nature of equation (jlj). The idea 
is to reduce the order of a higher order difference equation such that the analysis of the 
reduced system is feasible. Then, this is helpful when this analysis can facilitate the 
understanding of the dynamics for the original difference equation. 

Let fixi, X2, . . . , Xk+i) = '^^''+^ . Consider the vectorization of / i.e., 

Vf{xi, X2,..., Xk+l) = I — ^^^^ , Xi,...,Xk] . 
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Set H{xi, X2, ■ ■ ■ , Xk+i) = X1X2 . . . Xk+i and note that 



CXiX2---Xk+l cH(xi, X2, ■ ■ ■ , Xk+l) 

H{Vf{Xi,X2, . . . ,Xk+l)) - - 



1 + X1X2 . . . Xk+l 1 + H{xi, X2,..., Xk+l) ' 



So the map H makes Vf a (D, 0)-semiconjugate map with D being the nonnegative 
orthant of R.'''^^, i.e., D = [0, oo)'''^^, and the map 



i+t' 



serving as the factor map on [0, 00). Therefore, equation (jl]) is a semiconjugae factoriza- 
tion of the well-known Ricatti difference equation tn+i = (pitn) = In the sequence, 
we consider three cases as follow: 

Case /.■ < c < 1; In this case we claim that 

< Xn+i < c^"/''+^^ max{a;_fc, X-k+i, • • • , xo}, n > -k, 

where \n] is the greatest integer which is less than or equal to n. To prove the above 
claim note that it is true for all —k < n < —1. Assume that it is true for all integers 
less than or equal to some integer n{n > — 1). Then 

^Xfi—k 

Xn+1 = -7— 

i ~r XnXn—1 ■ ■ ■ Xn—k 

^ CXfi—k 

r n — fc -i I -1 

< c' '"^ max{x_fc,x_fc+i, . . . ,a;o| 
= c' ' max{x_fc,x_fc+i, . . . ,xo|, 

and this completes our inductive proof. Therefore, in this case every positive solution 
of equation (j4]) converges exponentially to origin. 

Case II: c > 1; In this case the factor map has a positive equilibrium t = c — 1. Re- 
call that in a Ricatti equation the positive equilibrium is globally asymptotically stable 
(see m. P79). Hence, t is globally asymptotically stable equilibrium of (p. Also, since the 
sequence {0(to)}J^O' ^0 £ (0, c>o) is bounded (this is evident by the global asymptotic 
stability of t), Theorem implies that the sequence {Vp{xo)}'^^Q, xq G (0,00)^"*"^ 

is bounded. Thereby, we conclude from Theorem 2(iv) that the sequence {Vp{xo)}'^^Q 
converges to the invariant fiber Hj^ since the sequence {0(to)}J^o converges to t. Now, 
consider the fiber Hj^. If (xi, X2, ■ ■ ■ Xk+i) G Hj^ then some simple algebra show that 

Vj~^^{xi, X2, . . . , Xk+l) = {xi, X2, ■ ■ ■ , Xk+l), 

Therefore, every member of Hj^ is a (A; + l)-cycle of V}; in other words, it is a periodic 
orbit of equation @ of period {k + 1). Then, in this case every solution of equation (jl]) 
converges to a (A; + l)-cycle. 
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Case III: c = 1; In this case origin is the unique equihbrium of equation (jlj) and (p. 
Some calculations show that 

0'(O) = 1, r(0) = -2, 

Therefore by Theorem 11.21 origin is semiasymptotic stable equilibrium of form the 
right, i.e., origin attracts the sequence {0"'(to)} for all tg G (0, oo). So, by an analysis 
precisely similar to that of used in the previous case we conclude, by Theorem 2(iv), 
that the the fiber Hq^ attracts every solution of equation @ since origin attracts every 
solution of the sequence {^"■(to)}- Now, assume that (xi, X2, ■ ■ ■ , Xk+i) € -f^o"^- Then 

X1X2 ■ --Xk+i = 0, 

and therefore, every solution of equation (jl]) either converges to origin or to a point in 
one of (A; + 1) coordinate planes in M^^^^, i.e., the following set 

S={{xi,X2,...,Xk+i)\Bl<i<k + l 3 x, = 0}. (7) 

elements of the attracting set 5* (note that S = USi, where Si are the coordinate planes) 
are (fc + l)-cycles of equation @. To prove this consider (xi, . . . , Xj, 0, Xj+i, . . . , Xk+i) G 
S. Then, it is easy to verify that 

Vf^^{xi, ...,Xi,0, Xi+i, Xk+i) = (xi, ...,Xi,0, Xi+i, Xk+i). (8) 

Consequently, every solution of equation (j4]) either converges to the origin or to a periodic 
orbit of period {k + 1), that is in a coordinate plane of M'^"'"^. We now summarize the 
above arguments into the following theorem that is one of the main result of this paper. 

Theorem 4.1. (Basic convergence Theorem) Assume Equation is given. Then, 

(i) For < c < 1, every positive orbit of equation ^ converges exponentially to origin. 

This implies that the equation has no other periodic orbit except the trivial fixed 
point (i.e., origin). 

(ii) When c > 1, any cycle of the equation is a {k + l)-cycle and the set of all 
{k + l)-cycles is given by equation Furthermore 

1. For = 1, every positive orbit of equation ^ either converges to origin or 
to a [k + l)-cycle. 

2. When c > 1, every positive solution of equation converges to a {k + 1)- 
cycle. 
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